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It is shown that the conductance of a superconductor-normal metal-superconductor junction 
can exhibit a significant dependence on the phase of the superconducting order parameter in the 
situation when the size of the normal region of the junction is much larger than the normal metal 
coherence length and the critical current of the junction is already exponentially small. The period 
of the conductance oscillations as a function of the phase can be equal to tt or 27r depending on 
parameters of the system. 

Suggested PACS index category; 05.20-y, 82.20-w 

^\ • The critical current of the superconductor-normal metal-superconductor (SNS) junction Ic ~ /coexp(— -^) decays 
exponentially and can be neglected when L >> Lt (See, for example, I^'). Here L is the length of the normal metal 
; I i: region of the junction shown in Fig. la, Lx = coherence length of normal metal, T is the temperature, 

' D = ^ is the electron diffusion coefficient, vp is the Fermi velocity and I is the electron elastic mean free path. On 
the other hand, the x dependent part of the conductance 6G{x) of the SNS junction can survive even in the case 
Lt <C L ^ Lin- Here x — Xi~X2, Xi,2 sue phases of the order parameters in superconductors composing the junction, 
', Lin = VDTin and Tin are inelastic diffusion length and inelastic mean free time respectively. The x-dependence of 
. SG originates from the fact that the amplitude of the Andreev reflection from the superconductor-normal metal (SN) 
boundary of an electron into a hole acquires an additional phase factor exp{ixi,2), while the amplitude of the reflection 
of a hole into an electron acquires the phase factor exp(— 1x1^2)- The weak localization contribution SGi{x) to '^^'(x) 
has been considered long ago . It arises in the first order approximation in the parameter ^ 1 and is connected 
with the interference of electrons traveling clockwise and counterclockwise along diffusive paths with close loops which 
contain Andreev reflections. Here pp is the Fermi momentum. The value of 6G1 is insensitive to the ratio between 
L and Lp and the characteristic energy interval which gives the main contribution into SGi is e ~ T. The period of 
SGi (x) as a function of x is tt [^1 . 

In this paper wc consider two other contributions SG2 and SG3 to 5G and show that the period of the oscillations 
of 6G as a function of x can be either tt or 27r depending on parameters of the system and the way how the 
conductance is measured. SG2 can be associated with the spatial coherence between electrons and holes arising due 
to Andreev reflection from the SN boundary ['^1. It arises in zeroth order approximation in the parameter ^ 1. 

^ . The contribution from this mechanism to the resistance of SN junction was considered in [^^i^l. It has been pointed 

O out that the clcctron-holc coherence is extended in the metal over the distance of order L^ — \ —■ Therefore, 
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it is clear that the main contribution to 5G2 comes from the relatively small energy interval e ^ Ec = ^ T, 
' ^t Lt ^ L <C Lin, SG2 decreases with L only as L~^. The period of (56*2 (x) as a function of x is 2tt. The 

, sum 6G1 + 5G2 gives the main contribution to 5G, provided the voltage drop V between the two superconductors 
^ composing the junction is zero and, consequently, x does not change in time. In this case, 5G(x) can be measured with 
\ the help of an additional lead "C" shown in Fig. la, while the phase difference x on the junction can be determined 
by an additional Josephson junction. In the case where the resistance of the SNS junction is measured by applying 
the voltage V between the superconductors, there is a third contribution 5Gz to 5G. The origin of SG^ is similar to 
the Debye relaxation mechanism of microwave absorption in dielectrics. In this case, due to the Josephson relation, 
X and, consequently, the quasiparticle density of states in the metal v{e, x) are functions of time. In other words, at 
small V the quasiparticle energy levels move slowly. The electron populations of the energy levels follow adiabatically 
the motion of the levels themselves and, as a result, the electron distribution becomes nonequilibrium. Relaxation of 
the nonequilibrium distribution due to inelastic processes leads to the entropy production, to the absorption of the 
energy of the external field and, therefore, contributes to 5G. 

We start with the calculation of 5G2- In the zeroth order approximation in the parameter the most adequate 
theoretical description of the system is provided in the framework of the Keldysh Green's function technique elaborated 
for superconductivity in [I2;i3,i4]^ ^-^c diffusive approximation and in the absence of electron-electron interaction 
in the normal metal region of the junction the linear response to the external electric field is described by the 
following set of equations: A). Usadel equations for the retarded normal g^{e,x) = cos9{e,x) and anomalous 
F^'{e,x) = — i exp(ix(e, a;)) sin0(e, a;) Green's functions 

^dle{t,x) + {ie-—)s\ne{e,x) ~ \(d^x{^,x)f sin2e{e,x) = 

Z Tin Z 
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describe the excitation spectrum of quasiparticles in the normal metal in contact with the superconductor. Here x 
is the coordinate along the junction (See Fig. la), and 0{e,x) = 9i{e,x) + i62{e,x) is a complex number. Inside the 
superconductor = f and 62 = 0. The boundary conditions for Eq.l have the form I^^l 

Dd,eie,x)=tcos{e{e,0+))cos{^ - ^^^) 

D sm{9ie, 0+))d^xie, x = 0+) = t sin(| - ^^^) (2) 

Here 0+ represents the normal metal side of the SN boundary, t = toVp and to is dimensionless transmission coefficient 
through the SN boundary. i3).The equation for the distribution function of quasiparticles /i, which describes the 
imbalance of populations of electron and hole branches of spectrum in metal has the diffusion form 

Dd^icosh" e2{e,x)dxfi{e,x)} = (3) 

The boundary conditions for Eq.3 have a form I^^l 

£>cosh02(e,O+)a,/i(e,O+) = i{/i(e,0+) - /i(e, 0")} sin(ei(e, 0+))cos-i(| - ^^^^) 

/i(e,a; = 0") = 0,/i(e,a; = L+) = -eVdJo{e) (4) 

Here /o(e) and V are Fermi distribution function and the voltage drop on the junction, respectively, 0^, L"*" stand for 
the superconducting sides of the SN boundaries. C). Finally the expression for the normal current densities across 
the junction is of the form; 

/+00 
cosh^ 62{e,x)dxh{e,x)de (5) 
-00 

Here vq is the density of states in the bulk normal metal. 

Using Eqs.(l)-(5) we get the expression for the resistance of SNS junction 

GsNS = GNLr dedJtan-^){ ^ + [ Jl dx'}'^ (6) 

i-oo ^ fcT^^cosh^2(e,0+)sin^i(e,0+)cos(|-^i^) io cosh^^aCe,^') ^ 

Here = <^dj^, ctd = e^Dv^ and S — Li x L2 are conductance of the normal metal part of the junction, Drude 
conductivity and the area of the junction respectively. The first and the second terms in Eq.(6) can be associated 
with the resistance of the SN boundary and the resistance of the normal region of the junction respectively. There 
are two major effects in metal, which arc due to proximity of the superconductor: 1) The effective diffusion coefficient 
in Eq.3 is renormalized due to Andreev reflection and is governed by the parameter 02. The correction to the 
local conductivity of the metal from this effect leads to the second term in Eq.(6). 2) The local density of states 
v{e,x) = v^Reg^^ {e.x) = i^o con 9i{e, x) cosh62{e, x) in the metal at small e is suppressed due to Andreev reflection 
and is governed by the parameter 9i. The contribution to the conductance of the SN boundary from this effect 
corresponds to the first term in Eq.(6). The x dependence of Gsns originates from the corresponding x dependence 
of 9\ and 92- Following to Eq.l, 2 near the SN boundary at small e the value of 0{e,x) should be close to its value in 
the superconductor 9s{e < A) = ^. It approaches its metallic value 9m = only after the distance L^. This means 
that the main contribution to Eq.(6) comes from the energy interval e ^ T. As a result, 

SG2 = aGN^dix) (7) 

Here g{x) is a universal function of x with the period 27r; a^la.tL^Lt~Y and a ^ ij^)^ at i <§; Lf. Nonlinerar 

effects of V are not important for the considered above mechanism, as long as, Ly = \J'^ <^ L 

We would like to mention that the interference corrections to the conductance, which are similar to Eq.7, has been 

measured [^^^1 in the Aharonov-Bohm type experiment. The geometry of the sample is shown in Fig. lb. In the case 
when L » Lt the amplitude of oscillations of the conductance of SN junction SGsn as a function of the Aharonov 
Bohm flux through the hole in Fig. lb is of order of 
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SGsN E, 27r$ 

Here is the flux independent part of the conductance of the sample, $o is the superconducting flux quanta and 
.92(2/) is a periodic function with the period 2-k. The temperature dependence of SGsn is in agreement with the 

experiment ^^^^ . 

Let us now discuss the Debye relaxation mechanism which arises in the case where the voltage V is applied between 
superconductors in Fig. la, and x changes in time 

? = ¥^ (9) 
dt h ^ ' 

Generaly speaking, in this case one has to solve a nonstationary version of Eqs.1,2. However, in the case when 
eV <C Ec one can use the adiabatic approximation where the time dependences of 0{€,x,xit)) and local density 
state i'{e,x,x{t)) originate from the corresponding time dependence of x{t)- The standard expression for the power 
absorption due to the Debye relaxation has the form (see, for example, I^^l) 

i'{e,x{t)) is the local density of states averaged over the volume v = LL1L2 of the normal metal region and breackets 
<> correspond to the averaging over the period of the oscillations Ti/eV . The main contribution to Eq.lO comes from 
the interval of energies e ^ Ec, where at x = the quasiparticle density of states is significantly suppressed compared 

with its bulk metal value. Using Eqs.1,2 at e <C E'c we get that P{e, x = 0) '^OyJ^ and D{e = 0, x{t)) ~ '^ol sin ^^|. 

When L <C Lt, the energy interval where i/(e) is small is of order of Et ~ indicating that the main contribution is 

coming from energy interval Et in this case. At Ec, Et ^ T one can neglect the e-dependence of Ti„{e) and, as a result, 
we have the expression for the contribution of this mechanism to the d.c. conductance of the junction [Q = V^SG^) 

SGs^aia^GN^^ (11) 

Li^ 

which can be even larger than Gjv. Here ai ~ 1 when Lt <C and ai ~ jfj^ when Lt » It is interesting that 
in the situation considered above the Debye absorption mechanism contributes to d.c. conductivity while in usual 
case it determines a.c. conductivity. Eq.(ll) is valid when ^^^p^ <C 1. In this limit one can introduce 5G:^{x{t)) 
which has the magnitude of order of Eq.ll and the period 2it. Following to Eq.lO, Q as a function of V saturates at 
eV^ which means that 6G3 decays with V as ( ^y^. )^- 

We would like to mention that the contribution of the Debye mechanism to the d.c. resistance of a close metallic 
sample of the Aharonov-Bohm geometry has been discussed in 1^*1. In that case the time dependence of the electron 
density of states was induced by the change of a magnetic flux $ through the ring. The important difference between 
the our case and the case considered in 1^*1 is that the average density of states in the normal metal is flux (and, 
consequently, time) independent. Therefore, the Debye absorption is nonzero only due to mesoscopic fluctuations of 
the density of states, whereas in the case of SNS junction, the; average density of states can be time dependent. Using 
results obtained in 1^*1 in the case when ^ 5o, Lt^ L, Lin ^ L,L2,L2 we can estimate the contribution to 5G 
due to the mesoscopic part of the Debye relaxation mechanism as 
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6GT - a2^Ec5orl (12) 



Here is the mean spacing between energy levels in the normal metal, a2 ~ 1 when Lt <C Lin and 02 ^ {^^}^ when 
Lt ^ Lin- We will neglect this contribution because, as we will see below, 5G^ <C SGi. 

Another contribution to 5G which arises in the first order approximation in the parameter is the above mentioned 
weak localization correction to the conductance of the junction SGi . The standard expression for the weak localizatin 
correction has the form SGi = —^^G^{f,f}, where C^{f,r') is, so called, Cooperon given by the diagram shown in 
Fig. 2, where solid lines correspond to electron Green functions, dashed lines correspond to impurity scatterings and 
triangles correspond to the Andreev reflections on the SN boundary. Since the Andreev reflections transform electrons 
into holes one has to introduce two kinds of Cooperons: electron and hole C'^ ones. We will use the following 
boundary conditions for G*^ and G'* , 
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t{C'{x,x') - eM'^iXi,2)C\x,x')) = Dd,C'{x,x')\,=o,L, 

d,C'ix,x') = d,C''{x,x'%=o,L (13) 

Eq.13 reflects the fact that in course of each Andreev reflection amplitudes of diffusion electron paths aquire mentioned 
above additional phases ±xi,2, but they do not take into account the spatial coherence between the electron and the 
hole which arises due to Andreev reflection. This is correct if L^, <^ L ov <^T. Here e* ^ T is the characteristic 
energy in the problem. As a result, the magnitude of the x-dependent part of the junction's conductance is of order 
of 

2 ( EcTin for OD case 

SGi = -a2^93{x) I ^ for ID case (14) 

[ \n{Lin/l)^ for 2D case 

Here g^{x) is a periodic function with the period x'- OD case corresponds to L, Li, L2 <C Li^', ID case corresponds to 
Lx ^ Lin ^ L2,L and 2D case corresponds to Li,L2 » I'm ^ L. The magnitudes of nonlinear in V effects due to 
the weak localization corrections depend of the dimensionality. In OD, ID, and 2D cases the requirement for linearity 

is eV<$:^. 

Tin 

Ratios between the three above considered contributions to 6G depend on the parameters and the dimentionality 
of the system. For example, in OD case at < eVrin/h we have 

6G1 ^ a2 ^_ . 5G^ ^ (I5) 
5G2 OL UGn SGs a\a^ TiGn Ec 



At large enough eV ^ ^ :p^Ecaia2 (but still smaller than Ec), SGs becomes much smaller than dG2- In this case, 

X dependent part of the resistance is determined by the sum {6G1 + 6G2)- For example, in OD case the ratio between 
6G1 (with the period tt) and 6G2 (which has the period 2n) is of the order of 

^-^^^ (16) 
6G2 a V UGn 

If y = and the conductance is measured with the help of the contact "C" in Fig. la, SG is the sum of 6G1 and SG2, 
ratio of which is determined by the correponding terms in Eq.(15),(16). 

These ratios can be both larger and smaller than unity, which means that the period of the oscillations of SG{x) 
can be either tt or 27r. This can explain why some experiments demonstrate tt periodicity of 5G I^'^'^o] t^]j\\q ^j^q others 
demonstrate 27r period I^^'^^l. The reason why 6G1, which arises only in the first order approximation in the small 
parameter ^ , can be comparable with 5G2 is that SG2 is determined by the small energy interval e <^ Ec <^T while 
5Gx is determined by e ~ T. 

All the three above considered contributions to 5G monotonically increase as T decreases. On the other hand, the 
conductance measured by four probe method can exhibit nonmonotonic behaviour as a function of T . The origin of 
this effect is that both in the bulk of the normal metal and in the vicinity of the SN boundary O2 is close to zero and 
the effective diffusion coefficient in Eq.3 is close to D, while at x ~ Lt we have ^2^1 and the effective diffusion 
constant can be a few times largrer than D. 

In the case V ^0 the current through the junction J{t) is described by Eq.5 and 

J = {G%j,s + mx))V{t) (17) 

Here G^gj^g is the x independent part of the conductance of the junction. Eq.l7 is valid when J 3> Ic- It follows from 
Eq.l7 that in the case when V{t) = Vq + V\ s,m.{p.t), the I-V characteristic of the junction averaged over the period of 
the oscillations has Shapiro steps at 4eVb = Ti^ '^1 or 2eVo = Ml depending on whether the main contribution to 5G 
comes from 5Gi or from 5G2 + SG3. The nonstationary effects discussed above can be observed, provided the frequency 
of the oscillations eV is much larger than the frequency of the phase slips in the junction uj ~ ^(|)^ exp(— ■^) 

which occur due to the thermal activation over the Josephson barrier. Here Ej = ^ is the Josephson energy of the 
junction. In conclusion, we would like to mention that the above presented results are valid if SgTin <C 1. 

We acknowledge usefuU discussions with B.Altshuler, H.Courtois, M.Devoret, D.Esteve, B.Pannetier, M.Sanquer 
and B.V.Wees. This work was supported by Division of Material Sciences, U.S. National Science Foundation under 
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